Edge states in graphene in magnetic fields 
— a speciality of the edge mode embedded in the n = Landau band 
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While usual edge states in the quantum Hall effect(QHE) reside between adjacent Landau levels, 
QHE in graphene has a peculiar edge mode at E — that reside right within the n = Landau level 
as protected by the chiral symmetry. We have theoretically studied the edge states to show that 
the E = edge mode, despite being embedded in the bulk Landau level, does give rise to a wave 
function whose charge is accumulated along zigzag edges. This property, totally outside continuum 
models, implies that the graphene QHE harbors edges distinct from ordinary QHE edges with their 
topological origin. In the charge accumulation the bulk states re-distribute their charge significantly, 
which may be called a topological compensation of charge density. The real space behavior obtained 
here should be observable in an STM imaging. 

PACS numbers: 68.37.Ef, 73.43.-f 
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Introduction — Ever since the anomalous quantum 
Hall effect (QHE) was experimentally observed, [2, HI fas- 
cination with graphene is mounting. The interests have 
been focused on the "massless Dirac" dispersions around 
Brillouin zone corners (K, K') in graphene, where the 
Dirac cone is topologically protected due to the chiral 
symmetry Q. The peculiar dispersion is responsible for 
the appearance of the n — Landau level (n: Lan- 
dau index) precisely around energy E = in magnetic 
fields. For the ordinary integer QHE an important gen- 
eral question is how the bulk and edge QHE conductions 
are related for finite samples. Many authors have ad- 
dressed this question, 0, Hj where one of the present au- 
thors has shown that the bulk QHE conductivity, a topo- 
logical quantity, coincides with the edge QHE conductiv- 
ity, itself another topological quantity. This constitutes 
a typical example of phenomena that, when a bulk sys- 
tem has a topological order d, 0,0, Q that reflects the 
geometrical phase of the system flOf. this should be re- 
flected and become visible in the edge states in a bounded 
system. [ll|, EH For graphene, two of the present authors 
and Fukui have shown that this "bulk-edge correspon- 
dence" persists in graphene, with both an analytic treat- 
ment of the topological numbers and numerical results 
for the honeycomb lattice [l3l[T3|. 

Now, in the physics of graphene, it is important to 
distinguish between the properties that arise from the 
continuum theory (i.e., the massless Dirac dispersion that 
comes from the k ■ p perturbation in the effective-mass 
formalism) from the properties that can only be captured 
by going back to the honeycomb lattice. In ref.[13l Il4j. 
we have already recognized this in a change from the 
Dirac to fermionic behaviors at van Hove singularities of 
the honeycomb lattice, and in the QHE edge modes that 
depend on whether the edge is zigzag of armchair. 

The purpose of the present paper is to reveal the 
features in the real-space profile of the edges states in 
graphene in magnetic fields B in the one-body problem. 
The graphene edge states in fact turn out to behave un- 



usually. A crucial point is that we find, from numeri- 
cally obtained dispersion and wave functions, that the 
E = edge mode, despite being embedded right within 
the n — bulk Landau level in the energy spectrum, has 
a wave function whose charge is accumulated along zigzag 
edges. This situation is drastically different from the or- 
dinary QHE where edge modes reside, in energy, between 
adjacent Landau levels, and their charge is depleted to- 
ward an edge. We can indeed realize in Fig. [1] (Landau 
spectrum vs position) that E = mode in graphene is 
special. The physics here points to a topological origin 
in a honeycomb lattice, which is in fact totally outside 
continuum models. In the absence of magnetic fields, a 
zigzag edge in graphene has been known to have a flat dis- 
persion at E = 0,[l6| which is protected by the bipartite 
symmetry of the honeycomb lattice. Here we are talking 
about the edge states in strong magnetic fields, which 
has a flat dispersion at E = 0. The charge accumulation 
along zigzag edges only occurs for the E = edge mode 
in the n = Landau level. To be more precise, the bulk 
states re-distribute their charge significantly on top of the 
zero-mode contribution, which may be called a topolog- 
ical compensation of charge density. On a practical side 
the present result predicts how an STM imaging should 
look like for graphene edges |15j. For comparison we have 
also examined edges in a bilayer graphene. [2fjl l2l| 



Single Layer — We consider the standard tight-binding 
model on the honeycomb lattice with nearest-neighbor 
hopping, where the magnetic field is introduced as a 
Peierls phase in the Landau gauge. The magnetic field 
is characterized by the flux in units of the magnetic 
flux quantum, <j) = BSq/(2tt) = 1/q in each hexagon 
with area 5*6 = (3v / 3/2)a 2 . Since we want to look at 
edges, we should be more explicit about the Hamilto- 
nian. Since honeycomb is a non-Bravais, bipartite lattice 
with two sublattice sites • and o per unit cell, we can 
define two fermion operators c. and c Q . For an armchair 




(a) Ordinary QHE 



(b) Graphene QHE 



FIG. 1: Schematic Landau-quantized spectra against real 
space position (x) for finite systems for ordinary QHE sys- 
tems(a) and graphene QHE(b). 




(a) Armchair edge 



(b) Zigzag edge 



FIG. 2: Honeycomb lattice with armchair (a) or zigzag (b) 
edges with ei , e2 being respective unit translation vectors. 
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FIG. 3: (Color online) Energy spectra against k?, (momentum 
along the edge) for a single-layer graphene in a magnetic field 
of = 1/5 with zigzag edges. Shaded regions are the bulk 
energy spectra, while red (blue) lines are the modes localized 
on the zigzag (bearded) edge. 
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(a) Armchair edge 
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(b) Zigzag edge 



edge (Fig.[^a)) the Hamiltonian then reads 
Ha =t ^2 [cl(j)c (j)+4(i + ei )c (j) 



c£(j +e a + e 2 )c.(j) 



H.c. 



Here j = jiei + J2&2 with ei,e2 defined in Fig. |U[a) 
specifies the position of a unit cell. For a zigzag edge 
(Fig. (Hb)) the Hamiltonian reads 

n * = 1 E [ c I(j)co(j) + e i;27rte ct(j) Co (j - ea) 



+4(j + ei)co(j) 



H.c. 



with ei,e2 as defined in Fig. Hfb). Hereafter we take t 
as a unit of energy and a as a unit of length. 

We assume that the system has left and right edges 
with a spacing L±, taken to be large enough [L\ = 5q 
here) to avoid interference. The length along the di- 
rection (e 2 ) parallel to the edge is also assumed to be 
long enough (L2), for which we apply the periodic bound- 
ary condition. We can then make a Fourier transform 

in that direction, c Q (j) = L^ 1 ^ 2 J2k 2 e%k2 ^ c a{hi ^2), for 
]2 = 1,2,- •• ,L 2 and a = »,o. This yields a k 2 - 
dependent series of one-dimensional Hamiltonian, Ti = 
y^ jk 7^id(/c2). The resultant eigenvalue problem reduces 



FIG. 4: Energy spectra against ki (momentum along the 
edge) for a single-layer graphene in a magnetic field of cf> = 
1/21 for armchair(a) or zigzag(b) edges. 



to HiD(k2)\ip(k2, E)) = E\tp(k2,E)), with corresponding 
eigenstates \t/j(k 2 ,E)). 

Having STM images in mind, we define the local charge 
density, 



2- 



Eo 



dE dk 2 \ME,Ji,k 2 )\ 2 . (1) 



Here x is the distance from the edge (as related to j\ via 
ei which is not normal to the edge), and E\ < E < E 2 is 
the energy window to be included in the charge density 
(which is normalized to unity when the window covers 
the whole spectrum). 

We have stressed that the E = edge mode is em- 
bedded within the n = Landau level, which is depicted 
in Fig. [31 a blowup of the energyspectrum (for a rela- 
tively high 0=1/5 for clarity). [13[ The shaded regions 
represent the n = Landau band, while the red curves 
the edge modes localized along the zigzag edge. We can 
see that, despite the presence of a strong magnetic field, 
there exists an exactly E = edge mode piercing the 
n = bulk Landau band. We can realize its topological 
origin by noting that there are an odd number (2q — 1) 
of edge modes with zig-zag edges, so that the bipartite 



(a) Armchair edge (b) Zigzag edge 



x/Ib x/Ib 
(a) Armchair edge (b) Zigzag edge 



FIG. 5: (Color online) Local charge density (oc area of 
each circle) for the single layer graphene with armchair(a) 
or zigzag(b) edges in a magnetic field <j> = I/2I for the energy 
window —0.05 < E < 0.05 around n = Landau level (see 
Figg). 



FIG. 6: (Color online) Scaled plot of the charge density I(x) 
against x/Ib, the distance from the edge normalized by the 
magnetic length, for the n = Landau level (marked with A in 
Fig- EJ with various values of magnetic field <j> for armchair(a) 
or zigzag(b) edges. 



symmetry (that forces an electron-hole symmetric energy 
spectrum) precisely dictates that the central edge mode 
has to be flat and at E = O.pl 

Figure^] shows the energy spectrum in a magnetic field 
<f> = 1/21 adopted hereafter for the armchair and zigzag 
edges. For this more moderate field the n = Lan- 
dau level around E = 0, with a narrow energy width 
(~ 0.05), almost looks like a line spectrum on this en- 
ergy resolution. We calculate the local charge density 
defined in eq.|T]) for the armchair and zigzag edges with 
the energy window —0.05 < E < 0.05 set to cover the 
n = Landau level (along with the embedded E = 
edge mode).[l9| In the result in Fig. [5] the charge density 
for an armchair edge decreases monotonically toward the 
edge, where the depletion occurs on the magnetic length 
scale (Ib — 3 3 / 4 a/\/27n/>)), as in ordinary QHE systems. 
In sharp contrast, a zigzag edge has the charge density 
for the •-sublattice that is accumulated toward the edge 
while the charge density for the o-sublattice is depleted. 

The question then is how the accumulation of the 
charge around the zigzag edge scale with the magnetic 
field. We plot in FigJS]the charge density I(x) normalized 
by the bulk value Jo (which is </>, when each Landau level 
of massless Dirac particles is fully occupied against 
the distance from the edge x measured by the magnetic 
length Ib for various values of the magnetic field with 
the energy window fixed at —0.05 < E < 0.05. In this 
scaled plot the profile of the zigzag edge states for vari- 
ous values of the magnetic fields fall upon common lines, 
where the accumulation of the charge on • sublattice as 
well as the depletion on o sublattice are seen to occur, re- 
spectively, on the magnetic length scale toward the edge. 
We have seen in Fig[3]that the flat dispersion for B / 
on a zigzag edge only exists over 1/3 of the Brillouin 
zone (that satisfies |1 - (-l) q e tqk2 \ < 1). We indeed ob- 
serve that the accumulated • charge over the depleted 
o charge, estimated by integrating the density over one 
half the sample width is (1 — <^>)/3 within the numerical 
accuracy. 

In order to confirm that the charge accumulation 
around zigzag edges is specific to the n = Landau level 
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FIG. 7: Scaled plot of the charge density I(x) against x/Ib 
for the n — 1 Landau level (marked with B in Fig. [4} with 
magnetic field (f> — 1/41 for armchair(a) or zigzag(b) edges. 



which embeds the edge mode, we can look at the charge 
density for n — 1 Landau level. The result in FigQhas 
I(x) monotonically decreases toward the edge for both 
armchair and zigzag edges, although we can notice that 
the charge density exhibits plateaus for each of • and o 
sublattices in a zigzag edge. 

Double Layer - We finally examine the bilayer 
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(a) crystal structure 
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(b) Energy spectrum (zigzag edge) 



FIG. 8: (Color online) (a) A bilayer graphene with Bernal 
stacking with the top layer having a zigzag edge. The transfer 
energies considered in the Slonczewski-Weiss-McClure model 
are displayed . (b) Energy spectrum for a bilayer graphene 
having a zigzag edge in the top layer with interlayer couplings 
t\ — t2 = O.lt in a magnetic field of <f> = 1/21. 
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(a) Armchair edge (b) Zigzag edge 



FIG. 9: (Color online) Local charge density (oc area of each 
circle) for a bilayer graphene with ti = t% = O.lt for arm- 
chair^) or zigzag(b) edges (indicated by arrows) in a mag- 
netic field 4> = 1/21 for the energy window —0.05 < E < 0.05. 
Red (blue) circles represent top (bottom) layer. 

graphene (Fig.[5[a)), which is interesting in its own right 
as many papers have pointed out, but the system is prac- 
tically interesting as well, since, experimentally, the STM 
imaging may be easier for the edge of the top layer resid- 
ing on a wider bottom layer. We consider the bilayer 
graphene with the AB (Bernal) stacking in the stan- 
dard Slonczewski-Weiss-McClure model[22|, HH, where 
there are two types of interlayer transfers, t\ and ti. 
For simplicity we have taken t\ = t% = O.lt, which are 
roughly the estimated values 0, HE, [26J. As in the sin- 



gle layer, the periodic boundary condition is applied in 
e2-direction, while a wider width (L\ = 7q) is taken for 
the bottom layer (with L\ = 5q for the top layer). De- 
spite the interlayer coupling, the Landau-quantized en- 
ergy spectrum (Fig. EUb)) is similar to those for the 
n = Landau level on this energy scale. Figure [H] dis- 
plays the charge density for armchair and zigzag edges 
with an energy window that covers the bilayer n = 
Landau level. We can see that the edge states in the 
top layer are similar to those in the single layer, namely, 
the charge density is accumulated toward the zigzag edge 
on one sublattice. So we can predict that a bright edge 
should be observed when an STM study is done for a 
zigzag edge of the top layer in a bilayer graphene. 

To summarize, we have shown that the charge density 
in graphene in strong magnetic fields should be totally 
unlike ordinary QHE systems, where the charge is accu- 
mulated toward zigzag edges. We wish to thank Hiroshi 
Fukuyama and Tomohiro Matsui for illuminating discus- 
sions, and for pointing out that bilayer graphene may 
be suitable for STM imaging. This work has been sup- 
ported in part by Grants-in-Aid for Scientific Research, 
No.20340098, 20654034 from JSPS and No. 220029004 
on Priority Areas from MEXT for MA, No.20340098, 
20654034 from JSPS and No. 220029004, 20046002 on 
Priority Areas from MEXT for YH, No.20340098 from 
JSPS for HA. 
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